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Abstract 

We present a space-time model of the collision of two homoge¬ 
neous, plane impulsive gravitational waves (each having a delta func¬ 
tion profile) propagating in a vacuum before collision and for which 
the post collision space-time has constant curvature. The profiles of 
the incoming waves are kS(u) and I5(v) where k,l are real constants 
and u = 0, v = 0 are intersecting null hyper surfaces. The cosmological 
constant A in the post collision region of the space-time is given by 
A = —6 k l. 
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1 Introduction 


Finding the space-time structure after the collision of gravitational and/or 
electromagnetic waves is a difficult problem in general relativity due to the 
non-linearity of the field equations. The problem is simplified by specialising 
to impulsive and/or shock waves which are plane and homogeneous and then 
exact solutions can be found, with the Khan-Penrose pQ and Bcll-Szekeres 
[2] solutions among the most famous. Up to now no solution where a cos¬ 
mological constant appears after the collision of two homogeneous, plane, 
impulsive gravitational waves has yet been found. This paper provides such 
an exact solution. 

Impulsive light-like signals, i.e. signals travelling with the speed of light 
and having a delta function profile, are idealised models of more realistic 
light-like signals having a profile with a finite width such as, for instance, 
a burst of light-like matter or of radiation. In some situations they may 
provide solvable models to describe their interactions and are sometimes 
used as classical models of quantum phenomena. In black hole physics one 
has the examples of mass inflation [3], the limiting curvature principle [0], 
Hawking radiation and quantum fluctuations |5] and internal structure of a 
Schwarzschild black hole In general relativity an impulsive light-like 

signal exists whenever the Riemann curvature tensor of the space-time man¬ 
ifold exhibits a delta function term with support on a null hypersurface, with 
the latter representing the space-time history of the signal and across which 
the first derivatives of the metric tensor are discontinuous. This signal can 
be a thin shell of light-like matter, or an impulsive gravitational wave or a 
mixture of both [8]. Recently the cosmological constant has received much 
attention in connection with a possible description of dark energy. Such ex¬ 
otic matter is described by a perfect fluid with an equation of state for which 
the sum of the energy density and isotropic pressure vanishes. In this paper 
we examine the possibility that the collision of two impulsive gravitational 
waves will produce such exotic matter by adopting a mathematical point of 
view, i.e. by solving Einstein’s field equations with appropriate boundary 
conditions. No attempt is made to propose a physical mechanism. 

For most of the known wave collision models in general relativity the same 
field equations apply before and after the collision. However the choice of 
field equations before collision does not determine the choice of held equations 
after collision. This freedom offers an opportunity to explore new and po¬ 
tentially interesting models. The well known space-time model of a head-on 
collision of two homogeneous, plane impulsive gravitational waves, travel¬ 
ling in a vacuum, involves the assumption that the post collision region of 
space-time is a vacuum space-time. With this assumption the post collision 
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region is described by the Khan-Penrose [T] solution of Einstein’s vacuum 
field equations (for a derivation see [9j). We demonstrate here that if the post 
collision region of space-time is assumed to be a solution of Einstein’s field 
equations with a cosmological constant then an exact solution of these field 
equations can be found satisfying the same conditions on the null hypersur¬ 
face boundaries of the post collision region as the Khan-Penrose solution. In 
addition the cosmological constant can be expressed simply in terms of two 
parameters which label each of the incoming waves. The post collision region 
is a space-time of constant curvature and is thus curvature singularity-free, 
in contrast to the Khan-Penrose model. The solution derived here is not 
an extension of the Khan-Penrose solution since it has the property that if 
the cosmological constant vanishes then at least one of the incoming waves 
vanishes. The post collision model presented here can be explained in terms 
of a redistribution of the energy in the incoming waves and this is described 
in some detail. 

In section 2 the incoming plane, impulsive gravitational waves propagat¬ 
ing through a vacuum are introduced, the collision problem is specified (as a 
light-like boundary value problem) and the solution of Einstein’s field equa¬ 
tions with a cosmological constant in the post collision region is given. This 
is followed in section 3 by a detailed study of the physical properties of the 
products of the collision which, in addition to a cosmological constant, in¬ 
clude impulsive gravitational waves (as in the Khan-Penrose collision) and 
light-like shells of matter. When reasonable physical restrictions are invoked 
the post collision region of space-time is anti- de Sitter space-time in this 
case. 


2 Colliding Waves 

A plane, homogeneous gravitational impulse wave propagating in a vacuum 
is described in general relativity by a space-time with line element 

ds 2 = —(1 + ku + ) 2 dx 2 — (1 — ku + ) 2 dy 2 + 2 du dv , (2.1) 

where k is a constant (introduced for convenience) and u + = ud{u ) where 
d(u) = 1 for u > 0 and d(u) = 0 for u < 0 is the Heaviside step function. The 
metric given via this line element satisfies Einstein’s vacuum field equations 
everywhere (in particular on u — 0). The only non-vanishing Newman- 
Penrose component of the Riemann curvature tensor on the tetrad given via 
the 1-forms d 1 = (1 + ku + )dx , d 2 = (1 — ku + )dy , d 3 — dv , d A = du is 

T 4 = -fc<5(u). (2.2) 
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Thus the curvature tensor is type N (the radiative type) in the Petrov clas¬ 
sification with the vector field d/dv the degenerate principal null direction 
and therefore the propagation direction of the history of the wave (the null 
hypersurface u — 0) in space-time. The wave profile is the delta function, 
singular on u — 0, and thus the wave is an impulsive wave. There are two 
families of intersecting null hypersurfaces u = constant and v = constant 
in the space-time with line element (12.ip . A homogeneous, plane impulsive 
gravitational wave propagating in a vacuum in the opposite direction to that 
with history u = 0 has history v = 0 and this is described by a space-time 
with line element 

ds 2 = —(1 + l v + ) 2 dx 2 — (1 — l v + ) 2 dy 2 + 2 du dv , (2.3) 

where l is a convenient constant and v + = vd(v). The Ricci tensor vanishes 
everywhere when calculated with the metric tensor given by this line ele¬ 
ment. The only non-vanishing Newnran-Penrose component of the Riemann 
curvature tensor on the tetrad given via the 1-forms tP = (1 + lv + )dx , d 2 = 
(1 — l v + )dy , id 3 = dv , $ 4 = du is 

= ~lS{v), (2.4) 

indicating a Petrov type N curvature tensor with degenerate principal null 
direction d/du. 

For the collision problem we envisage a pre-collision vacuum region of 
space-time v < 0 with line element (12.ip and a pre-collision vacuum region of 
space-time u < 0 with line element (12. 3 p (with both line elements coinciding 
when v < 0 and u < 0). The waves collide at u — v — 0 and the post 
collision region of the space-time corresponds to u > 0 and v > 0. In this 
region the line element has the form (in, ma, nu) 

ds 2 = —e~ u (e v dx 2 + e~ v dy 2 ) + 2 e~ M du dv , (2.5) 

where U, V, M are each functions of u, v. These functions must satisfy the 
following conditions on the null hypersurface boundaries of the region u > 
0 , v > 0: 

v = 0 , u > 0 =>- e~ u = 1 — k 2 u 2 , e v = - + ) M — 0 , (2.6) 

1 — k u 

and 

u = 0 , v > 0 e~ u = 1 — l 2 v 2 , e v = ^ \ V , M = 0 . (2.7) 

1 — l v 
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Einstein’s field equations with a cosmological constant A in the region u > 
0, v > 0 calculated with the metric tensor given by the line element (j2.5p 
read: 


Uuv 

= U U U V — A e~ M , 

(2.8) 

2 V 

^ v uv 

— U u V v + U v V u , 

(2.9) 

2 TJ 

n ^ uu 

= U 2 + V 2 -2M u U u , 

(2.10) 

2 U vv 

= U 2 + V 2 — 2 M v U v , 

(2,11) 

2 M uv 

= V u V v -U u U v , 

(2.12) 


where the subscripts denote partial derivatives. To implement our strategy 
below for solving (I2.8IWI2.12I) subject to the boundary conditions (I2.6j) and 
(12.71) we will need to know V v at v — 0, which we denote by (V v ) v= o, and V u 
at u = 0, which we denote by ( V u ) u= q . We already have from (12.61) and (12.71) : 


(Vv')v=o = - _ 2 ^2 and ( V v)u =o = j , (2-13) 

and also 

2 k 2 u 2l 2 v 

(U u ) v= o = ^ ^ 2 u 2 and \U v )u= o = i _ p v 2 ■ (2-14) 

In order to compute (V v ) v=0 and (V u ) u=0 we must first calculate (U v ) v=0 and 
(U u ) u= o- We obtain these latter quantities by evaluating (12.8|) at u = 0 and 
at v = 0 and solving the resulting first order ordinary differential equations. 
The constants of integration which arise are determined from the fact that 
U v and U u both vanish when u = 0 and v — 0, which follows from (12.141) . 
We then End that 


{Uv)v =0 


Aw (1 


|/c 2 m 2 ) 


k 2 u 2 


and (U u ) u= o = 


Au(l- \l 2 v 2 ) 
1 — l 2 v 2 


(2.15) 


Now evaluating (j2.9j) at v = 0 and at u = 0 provides us with a pair of first 
order ordinary differential equations for (I4)u=o an d (I4)u=o- These equations 
are straightforward to solve and the resulting constants of integration are 
determined from the fact that V u = 2 k and V v = 2 1 when u — 0 and v = 0, 
which follows from (12.131) . The final results are: 


iy^v')v=o 

(Tu) u=0 


21 + 


A 
3 k 


, A 
2k + Ti 


(i - ( tv )- 1 -' 2 - A 

V ; 3 k 


/1 + A; 2 w 2 \ 
\1 — k 2 u 2 ) 
1 + l 2 v 2 \ 

1 — l 2 v 2 ) 


(2.16) 

(2.17) 
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Dividing (j2.9j) successively by V u and by V v and then differentiating the 
resulting equations and combining them we obtain 


d 2 , V u _ ( TT K 
dudv ° S V„ V u V, 


uYa 

V V v/v 


(2.18) 


This suggests that we examine the possibility of a separation of variables: 


Vu _ A{u) 
K B(v) ’ 


(2.19) 


for some functions A{u) and B(y). The resulting mathematical simplification 
is that (12.19)) becomes a first order wave equation for V (see below) and that 
(12.18P becomes a second order wave equation for U. From a physical point 
of view we have shown |8] that if, as is the case in general, two systems of 
backscattered gravitational waves exist in the post collision region (one with 
propagation direction d/du in space-time and one with propagation direction 
d/dv) then (12.19)) implies that there exists a frame of reference in which the 
energy densities of the two systems of waves are equal. Using (12.13)) . (12.16)) 
and (12.17)) determines the right hand side of (12.19)) and the result is 


V u _ k [(l + VI - l 2 v 2 - ^(1 + l 2 v 2 ) 


Vv 1 K 1 + efl) VI - W - ^(1 + k 2 u 2 ) 

Hence this equation can be written as a first order wave equation 


( 2 . 20 ) 


V a = Vr, 


with u(u) and v(v) given by the differential equations 


du 

du 

dv 

dv 


1 + 


1 + 


A 

6 k l 
A 

6 k l 


Vl -k 2 U 2 --£-:(l + k 2 U 2 )] , 

6 kl 


( 2 . 21 ) 


( 2 . 22 ) 

(2.23) 


These two equations are interesting in general but we shall concentrate in 
this paper on two stand-out special cases: A = 0 and A = —6 kl. The 
case A = 0 is shown in the Appendix to correspond to the Khan-Penrose [T] 
space-time. 

With A = —6 kl we can solve (12.221) and (12.23)) . requiring u = 0 when 
u = 0 and v = 0 when v — 0, with 


u = tan 1 ku , v — tan 1 1 v . 


(2.24) 
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By (I2.2ip we have V — V(u + v) and the boundary condition (I2.6j) written 
in terms of u,v reads: when v — 0, V — log ( l_tan^ )- Hence 


V(u + v) 


log 


1 + tan (u + v ) 
1 — tan(h + v ) 


and restoring the coordinates u , v we have 


(2.25) 


V(u,v) 


log 


1 — k l u v + k u + l v' 
1 — kluv — ku — Iv 


(2.26) 


for u > 0, v > 0 provided A = —6 k l. Next writing (12. 9 p in terms of the 
variables u, v and using (12.211) and (12.251) we have 


Un + Uy 


8 tan(w + v) 

1 — tan 2 (h + v) 


which is easily integrated to yield 


e 


—U 


C(u 


/l — tan 2 (w + v)\ 
\1 + tan 2 (fZ + v) J 


(2.27) 


(2.28) 


where C(u — v) is a function of integration. When v = 0 the boundary 
condition (12.6p requires e~ u = 1 — tan 2 h and so C{u) = 1 + tan 2 u. Hence 
restoring the coordinates u,v we have U(u,v) given by 


u (1 — kluv) 2 — (ku + lv) 2 
(1 + k luv ) 2 


(2.29) 


for u > 0,n > 0. In the light of (I2.28p we see that U is a linear combination 
of a function of u — v and a function of u + v and thus satisfies the second 
order wave equation U U u = U m . This wave equation is the equation that 
(I2.18P reduces to when (12.211) holds and the barred coordinates are used. 

With V(u,v) and U(u,v ) given by (I2.26P and (I2.29P we use the field 
equation (12.81) with A = — 6kl to calculate M(u,v). The result is 


M (■ u, v ) = 2 log(l + k l uv ) , 


(2.30) 


and this clearly satisfies the boundary conditions (12.61) and (12.71) . Now with 
V, U and M determined a lengthy calculation verifies that the remaining field 
equations (12. 1 OH — (12 . 121) are automatically satisfied. Thus the line element 
(12.51) of the post collision region reads 


ds 2 = 


— (1 — kluv + ku + lv) 2 dx 2 — (1 — kluv — ku — lv) 2 dy 2 + 2 dudv 


(1 + kluv) 2 


(2.31) 
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If in the metric tensor components here we replace u, v by u + = u ^(u), v + = 
vd{y) we obtain in a single line element the expressions (12.ip and (12.3p for 
the pre-collision regions and (I2.3ip for the post collision region. In particular 
this will enable us to calculate the physical properties of the boundaries 
v — 0 , u > 0 and u — 0 , v > 0 of the post collision region. 


3 Post Collision Physical Properties 


With our sign conventions, choice of units for which c = G = 1, and energy- 
momentum-stress tensor T ab , Einstein’s field equations with a cosmological 
constant A read 

R ab = A g ab - 8 7t (T ab ~ -T c c g a b ) • (3.1) 

Also for a perfect fluid with proper density p and isotropic pressure p we 
have 

Tab = (p + p)u a u b -pg a b , (3.2) 

where u a , satisfying u a u a = 1, is the 4-velocity of a fluid particle. We 
note that for the exotic matter mentioned in section 1, p + p = 0 and thus 
T ab = p Qab and, with our sign conventions, Einstein’s held equations with this 
energy-momentum-stress tensor are equivalent to the held equations with a 
cosmological constant A = 8np. 

On the half null tetrad d 1 = e~^ u ~ v )d x ; p 2 — e~^ u + v )py ; $ 3 — 
e~^ M dv , d 4 = e~^ M du with V, U, M given by (12.261) . (I2.29j) and (12.301) with 
w, v replaced by u + , v + the Ricci tensor components of the space-time are 
given by 


Rab 


-6 k l d(u) i)(v) —— 6(v) SlSf 

2klv + (Pi’l~3) 4 4 

+ —Tnv— Hu)SA ' 


(3.3) 


This conhrms that the space-time region u > 0,v > 0 is a solution of the 
held equations with a cosmological constant, R ab = A.g a b, with A = —6 kl 
and that there are light-like shells with the boundaries v = 0,u > 0 and 
u — 0, v > 0 as histories, corresponding to the delta function terms in (13.31) . 
Here g ab are the (constant) metric tensor components on the half null tetrad 
given via the basis 1-forms {'$“}. The light-like shells have no isotropic 
surface pressure [12j and the surface energy densities are /qp and p( 2 ) given 
by 

An / k 2 u 2 — 3\ 

3 \1 — k 2 u 2 ) 
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8 vr /i(p 


on v = 0, u > 0 , 


(3.4) 















and 


(3.5) 


Ati /IV-3 

VK(2) = — 

The light-like shells must have positive surface energy densities. The only 
way to realise this on v — 0, u > 0 (respectively on u — 0,v > 0) is to 
have kl > 0 and k 2 u 2 < 1 (respectively kl > 0 and l 2 v 2 < 1). Thus the 
cosmological constant A = —6 k l must be negative. These restrictions on 
the coordinates are less restrictive than the condition k 2 u 2 + l 2 v 2 < 1 for 
u > 0 and v > 0 required in the Khan-Penrose post collision space-time 
on account of the presence of the curvature singularity. These restrictions 
on the coordinates also avoid infinite surface energy densities in the shells 
which are arguably as serious as a curvature singularity. Light-like shells did 
not appear in the Khan-Penrose model and their presence here is due to the 
non-zero cosmological constant. 

This result implies that the energy density of the exotic matter is nega¬ 
tive and that consequently only the so-called strong energy condition na, 
namely, p + p > 0 and p + 3p > 0, can be satisfied. 

The Newman-Penrose components of the Weyl conformal curvature ten¬ 
sor are given by 


on 


u — 0 , v > 0 


v bn = 


/(I + k 2 u\ 


■8(v) 


= 


k{ 1 + l 2 v\) 


5(u) , Ti = T 2 = T 3 = 0 


1 — k 2 u\ ’ 1 — l 2 v\ 

(3.6) 

Thus the boundaries v = 0, 0 < k 2 u 2 < 1 and u = 0, 0 < l 2 v 2 < 1 are the 
histories of impulsive gravitational waves corresponding to the delta function 
terms here. The post collision region u > 0, v > 0 is conformally flat and is a 
space-time of constant curvature with Riemann curvature tensor components 
given by 

Rabcd, 2 k l ip ad Qbc 9ac 9bd] • (3.7) 


Hence this region of space-time does not possess a curvature singularity, in 
striking contrast to the post collision region of the Khan-Penrose space-time. 


4 Summary 

We can briefly summarise our results as follows: for this model collision the 
energy in the incoming impulsive gravitational waves is re-distributed after 
the collision into two light-like shells of matter and two impulsive gravi¬ 
tational waves moving away from each other followed by a space-time of 
constant curvature. When the surface energy densities of the post collision 


9 





light-like shells of matter are required to be positive the space-time of con¬ 
stant curvature must be anti de Sitter space-time. 
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A The case A = 0 

With A = 0 we can solve (12.22)) and (12.23)) . requiring u = 0 when u — 0 and 
v = 0 when v — 0, with 

u = sm~ 1 ku, v = sin -1 1 v . (A.l) 
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By (I2.2ip we have V — V(u + v) and the boundary condition (I2.6j) written 

in terms of u, v reads: when v — 0, V — log ( |+ sm _ ). Hence 

0 V 1—sin it ) 


V(u + v) = log 


1 + sin (ri + v) 


= log 


1 — sin(w + v) / 
Restoring the coordinates u, v we have 

V (u, v ) = log 


cos u + sin v \ ( cos v + sin u 


cos u — sin v 


cos v — sin u 


\/l — k 2 u 2 + l v \ (y/1 — l 2 v 2 + k 


u 


y/1 — k 2 U 2 — l V ) \ y/1 — l 2 V 2 — k 


u. 


(A.2) 

(A.3) 


for u > 0, v > 0 provided A = 0. Writing (1 2.9 j) in terms of the variables u, v 
and using (I2.21j) and ()A.2jl we have 


Uu + Uy — 2 tan(h + v) . (A.4) 

Integrating this results in 

e~ l = D(u — v) cos (u + v) , (A.5) 

where D(u — v) is a function of integration. With v = 0 the boundary 
condition (12.61) requires e~ u = cos 2 u and so D(u) = cos u. Restoring the 
coordinates u,v we have U(u,v ) given by 

e~ u = cos (u — v) cos(h + v) = 1 — k 2 u 2 — l 2 v 2 . (A.6) 


Now (12.8p with A = 0 is automatically satisfied. Combining (12.8p with A = 0 
and (12.121) we have 

(2 M + U) uv — V u V v . (A.7) 

Changing the independent variables u, v to u, v using flA.ip . and using (IA.2p . 
this reads 

(2 M + U)uv = 4 sec 2 (h + v) . (A.8) 

Integrating and using (IA.6P we arrive at 


_ 2 m cos 3 (u + v) . 

e =- t - 7-f(u)g(v) , 


cos (u — V ) 


(A.9) 


where f(u) and g(v) are functions of integration. From the boundary condi¬ 
tions (12.6 p and (12.7p we see that M = 0 when 7 = 0 and M = 0 when 7 = 0 
and so it follows that f(u)g(v ) = sec 2 u sec 2 v. Hence 


e 


-M 


cos 3 (u + v) 

a -—-— sec u sec v . 

\ cos (u — v) 


(A.10) 
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Restoring the coordinates u, v and simplifying this becomes 


11 ) 

for u > 0,u > 0 when A = 0. Substituting (1A.3I) . (IA.6|) and flA.111) into 
the field equations (12.10[) and (ETTTj) verihes that these latter equations are 
automatically satisfied. The line element (12. 5 j) with V, U and M given by 
(I A. 31) . (1A.6|) and (1A.11I) is the Khan-Penrose post collision line element. No 
derivation is given in [I], 


e~ M = 


(1 — k 2 u 2 — l 2 v 2 ) 3 / 2 


(\/l — k 2 u 2 \/ 1 — l 2 v 2 + klu v) 2 y/l — k 2 u 2 \J 1 — l 2 v 2 


(A. 
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